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1 Abstract 

Effect of the trimodal random magnetic field distribution on the phase diagrams of the anisotropic 
quantum Heisenberg model has been investigated for three dimensional lattices with effective field 
theory (EFT) for a two spin cluster. Variation of the phase diagrams with the random magnetic 
field distribution parameters has been obtained and the effect of the anisotropy in the exchange 
interaction on the phase diagrams has been investigated in detail. Particular attention has been 
devoted on the behavior of the tricritical points with random magnetic field distribution. Keywords: 
Quantum anisotropic Heisenberg model; random magnetic field; trimodal distribution 

2 Introduction 

Recently there has been growing theoretical interest in the random field lattice spin models. For 

• instance, Ising model in a quenched random field (RFIM) has been studied over three decades. 
The model was introduced for the first time by Larkin [T] for superconductors and later general- 
ized by Imry and Ma [2]. Diluted antiferromagnets (such as Fe x Zn\- x Fi, RbiCo x Mg\- x F^ and 

, Co x Zn\- x F-2) in a homogenous magnetic field behave like ferromagnetic systems in the presence of 

• random fields [3J H] - Beside this, a rich class of experimentally accessible disordered systems can be 
described by RFIM, such as structural phase transitions in random alloys, commensurate chargc- 
density-wave systems with impurity pinning, binary fluid mixtures in random porous media, and 

\ the melting of intercalates in layered compounds, such as TiS^^- RFIM generally mimics the 
phase transitions and interfaces in random media [6] [7], e.g prewctting transition on a disordered 
substrate can be mapped onto a 2D RFIM problem [5]. Also, RFIM has been applied in order to 
describe critical surface behavior of amorphous semi- infinite systems [TU] . 

Random field distribution of the magnetic field produces drastic effects on the phase diagrams 
and related magnetic properties of the system. It has been shown that Ising systems under the 
influence of discrete symmetric distributions, like bimodal and trimodal |12| distributions, 
show tricritical behavior, while continuous symmetric distributions like Gaussian distribution [13) 
exhibit only second order transitions. 

On the other hand, as far as we know, there have been less attention paid on the random field 
effects on the Heisenberg model, which is more realistic model than the Ising model for the spin 
systems. Albuquerque and Arruda |14| studied the effect of the bimodal random field distribution 

\ on the phase transition characteristics of the spin-1/2 isotropic classical Heisenberg model and they 
found tricritical behavior within the EFT formulation for the two spin cluster which is abbreviated 

\ as EFT-2. Oubelkacem et al., studied the same system with another approach, namely EFT 
with probability distribution technique and they obtained similar results [15]. Albuquerque et 
al. [16] treated the same system with amorphisation effect, again with the EFT-2 formulation. 
Recently, Sousa et al. have studied the effect of the bimodal random field distribution on phase 
transition characteristics of the isotropic -classical and quantum- spin-1/2 Heisenberg model within 
the EFT-2 formulation and also they found a tricritical behavior [T7]. All these works have been 
restricted to the spin-1/2 isotropic Heisenberg model with bimodal random field distribution and 
they concluded that tricritical behavior exists in this system as in Ising model with bimodal random 
field distribution. They utilized an EFT which is characterized by differential operator technique 
introduced by Honmura and Kancyoshi for Ising systems |18j . EFT approximation can provide 
results that are superior to those obtained within the traditional mean field approximation, due 
to the consideration of self spin correlations which are omitted in the mean field approximation. 
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EFT for a typical Ising system starts by constructing a finite cluster of spins which represents the 
system. Callen-Suzuki spin identities QIJ1 HD] are the starting point of the EFT for the one spin 
clusters. If one expands these identities with differential operator technique, multi spin correlations 
appear, and in order to avoid from the mathematical difficulties, these multi spin correlations are 
often neglected by using decoupling approximation |21j . Working with larger finite clusters will 
give more accurate results. Callen-Suzuki identities have been generalized to two spin clusters in 
Ref. [23] (EFT-2 formulation). This EFT-2 formulation has been successfully applied to a variety 
of systems, such as quantum spin-1/2 Hcisenberg fcrromagnet |231 124 ] and antiferromagnet [25) 
systems, classical n- vector model [26] [27], and spin-1 Heisenberg ferromagnet [28l [29] . 

The aim of this work is to investigate the effect of the symmetric discrete random field distri- 
butions (bimodal and trimodal) on the phase transition characteristics of a spin-1/2 anisotropic 
quantum Heisenberg model on simple cubic and body centered cubic lattices. Quantum Heisenberg 
model can take into account the quantum fluctuations which dominates the thermal fluctuations in 
the low temperatures. Thus it is expected that it gives more reasonable results than the classical 
one at this low temperature region. We follow the EFT-2 formulation which is derived in Ref. [23] 
for this system. 

The paper is organized as follows: In Sec. [3] we briefly present the model and formulation. 
The results and discussions are presented in Sec. 0J and finally Sec. [5] contains our conclusions. 



3 Model and Formulation 

We consider a lattice which consists of TV identical spins (spin-1/2) such that each of the spins has 
z nearest neighbors. The Hamiltonian of the system is given by 

« = - E (^ s ? + + - E /7 << w 

<i,j> i 

where sf, s\ and sf denote the Pauli spin operators at a site i. J Xl J y and J z stand for the anisotropy 
in the exchange interactions between the nearest neighbor spins and Hi is the longitudinal magnetic 
field at a site i. The first sum is carried over the nearest neighbors of the lattice, while the second 
one is over all the lattice sites. Magnetic field is distributed on the lattice sites according to a 
trimodal distribution function which is given by 

P (Hi) = P 5 [Hi) + ^ [5 (Hi -H ) + 5 (H t + H )} (2) 

where p is a real number which provides < p < 1, and 6 stands for the delta function. The 
distribution given by Eq. ([2]) covers a bimodal distribution for p = and reduces to the system 
with zero magnetic field (pure system) for p = 1. According to the distribution given in Eq. ([2]), p 
percentage of the lattice sites are subjected to a magnetic field Hi = 0, while half of the remaining 
sites are under the influence of a field Hi = Hq whereas the field Hi = —Hq acts on the remaining 
sites. 

We use the two spin cluster approximation as an EFT formulation namely EFT-2 formulation[23j 
In this approximation, we choose two spins (namely si and S2) and treat interactions exactly in 
this two spin cluster. In order to avoid some mathematical difficulties we replace the perimeter 
spins of the two spin cluster by Ising spins (axial approximation) [24] . After all, by using the 
differential operator technique and decoupling approximation (DA) |21j . we get an expression for 
the magnetization per spin as 

m = /i (sf + s^ = ([A x + mB x ] Z0 [A y + mB y } z ° [A xy + mB xy ] Zl ) F (x,y, H ) \ x =o, y =o (3) 

where each of si and S2 has number of zq distinct nearest neighbors and both of them have z\ 
common nearest neighbors. The coefficients are defined by 

A x = cosh (J Z V X ) B x = sinh (J Z V X ) 

A y = cosh(J z V !/ ) B y = sinh(J z V y ) (4) 

A X y = COSh [J Z (V X + S7y)} B X y = SUlIl [J Z (V X + Vj,)] 



2 



where \7 X = d/dx and V a = d/dy are the usual differential operators in the differential operator 
technique. Differential operators act on an arbitrary function via 

cxp (aV, + b\7 y )G (x, y) = G(x + a,y + b) (5) 

with any constant a and b. The function in Eq. ((3]) is given by 

F{x,v,H )= fdH 1 dH 2 P(H 1 )P(H 2 )f(x > y,H 1 ,H 2 ) (6) 



where 

f( R R] = x + y + H 1 +H 2 sinh (PX ) 

J{X,y ' U 2> X Q cosh (/3X ) + exp (-2/3 J z ) cosh ((3Y ) [ ) 



and 



1/2 



(J x — J y ) + (x + y + Hi + H 2 ) , Y = (J x + J y y + ( x -y + H 1 -H 2 y . (8) 



1/2 



Here /? = l/(fcsT) where ks is Boltzmann constant and T is the temperature. 
With the help of the Binomial expansion, Eq. ([3]) can be written as 



ZQ Zq 21 

m = EEE c ^ mP+9+r ( g ) 

p=0(}=0r=0 

where the coefficients are 

C' pqr =( Z °^ Z °^j( Zl r ) ^-"^-^-S^'B:/^.!,.^) U=o,„=o (10) 

and these coefficients can be calculated by using the dchnitions given in Eqs. (j4]) and ([5]). Let us 
write Eq. (j9]) in more familiar form as 



= ]TC fe m fc (11) 



m 

Zq Zo Z\ 

C k = y^y^y^^ P+q+r,kC' pqr (12) 

p=0g=0r=0 

where 5ij is the Kronecker delta. It can be shown from the symmetry properties of the function 
defined in Eq. ((6]) and operators defined by Eq. (j4]) that for even k, the coefficient Ck is equal to 
zero. 

For a given set of Hamiltonian parameters ( J x , J y , J z ), temperature and field distribution pa- 
rameters (p,Hq), we can determine the coefficients from Eq. (|12p and we can obtain a non linear 
equation from Eq. (fTTj) . By solving this equation, we can get the magnetization (m) for a given 
set of parameters and temperature. Since the magnetization is close to zero in the vicinity of the 
critical point, we can obtain a linear equation by linearizing the equation given in Eq. {TTJ which 
allows us to determine the critical temperature. Since we have not calculated the free energy in 
this approximation, we can locate only second order transitions from the condition given as 

Ci = 1, C 3 < (13) 

The tricritical point at which second and first order transition lines meet can be determined from 
the condition 

d = 1, C* 3 = 0. (14) 
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4 Results and Discussion 

Let us scale the exchange interaction components with the unit of energy J as, 

r„ = — 
J 

where n = x,y, z. Let us choose r z = 1, then r x , r y can be used as the measure of the anisotropy 
in the exchange interaction. It can be seen from the definition of the function given in Eq. ([7]) 
that the transformation J x — > J y ,J y Jx does not change the function. Hence, we can say 
that mentioned transformation does not affect the formulation. Because of that let us fix r x and 
concentrate only on varying r y values. Our investigation is on the simple cubic (zq = 5, Z\ = 0) 
and body centered cubic (zq = 7, Z\ = 0) lattices. 



4.1 Bimodal Distribution 

This distribution is given by Eq. ([2]) with p = and it distributes the longitudinal magnetic fields 
±i?o to the lattice sites with equal percentages. Increasing randomness (which comes from increas- 
ing Hq values in bimodal distribution) will reduce the critical temperature, as expected. Beside 
this, in order to concentrate on the variation of the tricritical point as a function of anisotropy in 
the exchange interaction, let us examine the evolution of the phase diagrams in a (fcsT c / J, Hq/ J) 
plane with different r y values. In Fig. (fT]), we can see the phase diagrams of the anisotropic 
quantum Hcisenbcrg model on simple cubic and body centered cubic lattices for a bimodal ran- 
dom magnetic field distribution in a (ksT c / J, Hq/ J) plane for some selected values of the r y . We 
can track the path of the tricritical point (denoted by filled circles) with changing r y values for 
simple cubic and body centered cubic lattice in Figs. (fT]) (a) and (JlJ (b), respectively. Increasing 
r y values -which means that the anisotropy in the spin-spin interaction increases- reduce the two 
coordinates of the tricritical point in (ksT c / J 1 Hq/ J) plane. Also for arbitrarily fixed Hq/ J values, 
rising anisotropy reduces the critical temperature, as expected. For a certain r y value, we can see 
that tricritical point of the body centered lattice appears to be greater than that of the simple 
cubic lattice. We can also see from Fig. (fT]) that falling rate of a critical temperature corresponding 
to a certain Hq/ J value rises with increasing r y . Tricritical point coordinates of the isotropic case 
(r x = f y = 1) found in the present work (Hq/ J,kBT c / J) = (2.274,2.748) can be compared with 
the classical case (Hq/ J.ksTc/ J) = (2.389,2.785) [THOU]. Our values in the quantum case are 
slightly lower than corresponding classical values, as expected. 




Figure 1: Phase diagrams of the anisotropic quantum Hcisenbcrg model for (a) simple cubic 
and (b) body centered cubic lattices with the bimodal random magnetic field distribution in 
(ksTc/ J, Hq/ J) plane for some selected values of r y . Solid lines represent the second order tran- 
sitions, while the filled circles denote the tricritical points. Fixed parameter value is r x = 1.0. 
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4.2 Trimodal Distribution 

As we can see from Eq. ([2]) that trimodal distribution brings the system closer to the pure one 
(system with zero magnetic field) with increasing p. Thus, it is expected that the phase diagrams 
in {UbTc/ J, H / J) plane with increasing p will become parallel lines with the Hq/J axis (at the 
value of p = 1), which is nothing but just the critical temperature of the pure system with given 
parameter values. Also for a given Hq/ J value, increasing p means decreasing randomness effect. 
Therefore, the ferromagnetic region is expected to get expanded in {ksT c / J 1 Hq/ J) plane. However, 
it would be interesting to inspect again the evolution of the tricritical point with increasing p 
values. In Fig. (|2[). we can see the phase diagrams of the anisotropic quantum Heisenberg model 
in (kBT c /J,Ho/J) plane for simple cubic and body centered cubic lattices in the presence of 
a trimodal random magnetic field distribution, and for some selected values of p. We can see 
from Fig. that observed tricritical points at low p values gradually decrease and reduce to 
ksTc/J = at a certain Hq/ J. Right after the disappearance of the tricritical point, the system 
exhibits a second order reentrant behavior, i.e. after a certain value of p (which depends on the 
lattice geometry, as well as r x ,r y values) the system which stays in a disordered phase at zero 
temperature can pass to an ordered phase with a second order transition due to increasing thermal 
fluctuations, then it passes again to another disordered phase characterized by a second order 
transition with increasing temperature. For p values greater than a special value (let us denote it 
by p*), the phase diagrams in (fcsT c / J, Hq/ J) plane do not intersect the x-axis, i.e. for p > p* , the 
system will stay in an ordered phase at zero temperature for any Hq / J value. This specific value 
depends on the lattice geometry and the degree of the anisotropy in the exchange interaction (r x 
and r y ). When we compare Figs, ©(a) and (b) with each other, we can conclude that increasing 
exchange anisotropy does not alter this situation significantly for small anisotropy, but it affects 
only the critical values (e.g. compare Figs. ([2]) (a) and (b) or Figs. ([2]) (c) and (d)). The same 
observation also holds for a body centered cubic lattice. After all, we can say that in general, three 
qualitatively different regions arc observed in the phase diagrams depicted in a (fcsT c / J, Hq/ J) 
plane as p rises. Namely, a tricritical behavior region, a second order reentrant region, as well as a 
region which gets expanded towards increasing Hq / J direction. We also note that the second order 
reentrant region corresponding to the simple cubic lattice is wider than that of the body centered 
cubic lattice. This point is depicted more apparently in Fig. ([3]). 

In order to see the dependence of p* on the exchange anisotropy, we plot the variation of p* 
with r y for some selected values of r x in Figs. ([3]) (a) and (b) for simple cubic and body centered 
cubic lattices, respectively. We note that for a certain value of r x , the points which lie below the 
related curve in (p*, r y ) plane have disordered ground states for high Hq/ J values. In other words, 
for a selected r x value, if we choose (p, r y ) pairs which lie under the curve corresponding to selected 
r x , this means that the phase diagram of the system in (fcgT c / J, Hq/ J) plane either exhibits a 
tricritical behavior or a second order reentrant behavior. For a system with parameter values of 
(p, r v ) which lie above the related curve has an ordered ground state in high Hq/ J region, i.e. the 
phase diagram of the system in a (fcsT c / J, Hq/ J) plane does not intersect the x-axis. Thus, we can 
see from Fig. ([3]) that increasing anisotropy contracts the region where the phase diagrams on the 
right side of the (ksT c / J, Hq/J) plane are stretched. For a fixed value of r x , and after a certain 
value of r y , we get p* = 1.0 which means that the system can not exhibit an ordered ground state 
at high values of Hq/J regardless the value of p. Another point is that the simple cubic lattice has 
a higher p* value for a fixed r x ,r y , unless we have p* = 1.0. 

Finally, in order to reveal the relation between the p* value and the anisotropy in the exchange 
interaction more clearly, we depict equally valued p* curves in {r x ,r y ) plane in Fig. |4]). As we 
can see from Fig. |4]) that these curves are symmetric about the r x = r y line. This is obvious, 
since r x — > r y ,r y — > r x transformation does not alter the equations of the system, as stated in 
the beginning of this section. In general, for a selected value of r x , if we move on the increasing 
direction of r y , we pass through increasing p* curves, successively. But for instance, if we choose 
r x = 2.0 and continue in increasing direction of r y starting from r y = 1.0 (in Fig. ((4]) (a)), 
we intersect p* = 0.54 curve two times and the other curves once. Thus, the relation between 
anisotropy and p* is not linear; i.e. we can not say p* increases monotonically with increasing 
anisotropy. 
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Figure 2: Phase diagrams of the anisotropic quantum Heisenberg model for (a),(b) simple cubic 
and (c),(d) body centered cubic lattices with a trimodal random magnetic field distribution in 
(ksT c / J, Ho/ J) plane for selected values of p from p = 0.0 to p = 1.0 with increment of 0.1. 
Increment direction is shown as arrows in each figure. Solid lines represent the second order 
transitions, filled circles denote the tricritical points. Fixed parameter value is r x = 1.0. 

5 Conclusion 

The effect of the trimodal random magnetic field distribution on the phase diagrams of the 
anisotropic quantum Heisenberg model has been investigated in detail. The effects of the random 
magnetic fields in the presence of exchange interaction anisotropy have been discussed for three 
dimensional lattices, namely for simple cubic and body centered cubic lattices. Qualitatively sim- 
ilar characteristics have been observed for the phase diagrams of these different lattice geometries. 
However, quantitative differences have been obtained for bimodal and trimodal distributions. 

Calculations for bimodal distribution performed in this work may be thought as a generalization 
of the earlier works which were devoted to a bimodal distribution in isotropic Heisenberg model [1^1 
[Tbilisi E] case to the anisotropic case. According to the phase diagrams plotted in {k B T c j J, H /J) 
plane, presence of exchange anisotropy in the system reduces all critical temperatures, as well as 
tricritical points. 

On the other hand, trimodal distribution gives rise to qualitatively three different phase diagram 
types with varying p values in (ksT c / J, Hq/ J) plane. Namely, as p increases starting from p = 
0.0, we observe phase diagrams exhibiting tricritical behavior, phase diagrams with second order 
reentrant phenomenon, and phase diagrams which become stretched to the right hand side of 
(fcsT c / J, Hq/ J) plane, respectively At the end of this evolution process (i.e. for p = 1.0), phase 
diagrams become parallel lines with the Hq/J axis. In this case, fcsT c /J value just corresponds 
to the critical temperature of the pure system. This classification scheme is valid for both simple 
cubic and body centered cubic lattices. 

Particular attention has been devoted on the p value (namely p* ) for which the phase diagrams 
become stretched to the right hand side of {ksTc/ </, Hq/ J) plane. The effect of the anisotropy 
in the exchange interaction on p* value has been investigated in detail. For this purpose, equally 



G 




Figure 3: Variation of p* with r y for some selected values of r x for the anisotropic quantum 
Heisenberg model with a trimodal random magnetic field for (a) simple cubic, and (b) body 
centered cubic lattice. 




Figure 4: Equally valued p* curves in (r x , r y ) plane for the anisotropic quantum Heisenberg model 
with a trimodal random magnetic field on (a) simple cubic, (b) body centered cubic lattice. 



valued p* curves have been obtained in (r x ,r y ) plane. It has been shown that, the effect of the 
parameters r x , r y on the value of p* is not linear. 

We hope that the results obtained in this work may be beneficial form both theoretical and 
experimental point of view. 
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